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Abstract 

In a previous Letter, we outlined an approach to the calculation of quantum amplitudes 
appropriate for studying the black-hole radiation which follows gravitational collapse. This 
formulation must be different from the familiar one (which is normally carried out by 
considering Bogoliubov transformations), since it yields quantum amplitudes relating to the 
final state, and not just the usual probabilities for outcomes at a late time and large radius. 
Our approach simply follows Feynman's +ie prescription. Suppose that, in specifying the 
quantum amplitude to be calculated, initial data for Einstein gravity and (say) a massless 
scalar field are specified on an asymptotically-flat space-like hypersurface E/ , and final 
data similarly specified on a hypersurface E p , where both E j and E p are diffeomorphic 
to M"^ . Denote by T the (real) Lorentzian proper-time interval between Ej and E^? , as 
measured at spatial infinity. Then rotate: T — \T\ exp{—i9) , < 9 < 7r/2 . The classical 
boundary-value problem is then expected to become well-posed on a region of topology 
/ X M'^ , where / is the interval [O, |T|] . For a locally- supersymmetric theory, the quantum 
amplitude is expected to be dominated by the semi-classical expression exp(i5'ciass)5 where 
'S'ciass is the classical action. Hence, one can find the Lorentzian quantum amplitude from 
consideration of the limit 9 ^ Oj^ . In the usual approach, the only possible such final 
surfaces "Ep are in the strong- field region shortly before the curvature singularity; that is, 
one cannot have a Bogoliubov transformation to a smooth surface 'after the singularity'. 
In our complex approach, however, one can put arbitrary smooth gravitational data on 
E p , provided that it has the correct mass M ; thus we do have Bogoliubov transformations 
to surfaces 'after the singularity in the Lorentzian-signature geometry' — the singularity 
is simply by-passed in the analytic continuation (see below). In this Letter, we consider 
Bogoliubov transformations in our approach, and their possible relation to the probability 
distribution and density matrix in the traditional approach. In particular, we find that 
our probability distribution for configurations of the final scalar field cannot be expressed 
in terms of the diagonal elements of some density-matrix distribution. 

1. Introduction 

In [1-4] , based on the thesis [5] , we described a calculation of the quantum amplitude 
for an initial spherically-symmetric configuration of gravity and a massless scalar field (p 
to make a transition to another nearly-spherically-symmetric configuration at a very late 
Lorentzian time T, apart from some weak perturbations ih\^\ ^*^^'')f • The boundary data 
are specified on an initial and a final asymptotically-fiat hypersurface E/ , and the time 
T is measured at spatial infinity. The Lorentzian-signature space-time metric g^i, {iJ,,u = 
0, 1, 2, 3) has spatial components hij = gij (i.j = 1, 2, 3) on surfaces with t = x'^ = const.. 
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such as S/ and Hp • This transition could refer to the quantum evaporation of a black hole, 
following gravitational collapse from a slowly-moving, diffuse, spherical configuration. The 
final configuration could refer to the nearly-spherical stream of (massless) radiation, and its 
corresponding Vaidya-like gravitational field [6,7], as viewed on a space-like hypersurface 
Ei? , of topology IR^, which cuts through all the radiation. The weak perturbations on Ei? , 
in the present field language, correspond to emitted particles. 

The calculation proceeds following Feynman's +ie prescription [8]. If T is rotated 
into the complex: T — > |T|exp(— z^), for < ^ < 7r/2, then the c/assica/ boundary- 
value problem, given the above nearly-spherically-symmetric boundary data, is expected 
to become well-posed [1]. 

For a simple example (cf. [10]), consider the 2-dimensional Euclidean-space Laplace 
equation 

^ + ^=0; 0<T<T, -oo<a;<+oo, (1.1) 

subject to boundary data of the form (say) 

(/.(r = 0, a;) = , (/.(r = T, x) = M^) , (1-2) 

on the 'initial and final surfaces'. Here, the boundary data given by 4>i{x) should be 
thought of as a C°° function, of rapid decrease as |x| — > oo . 

Let $i(/c) denote the Fourier transform of (j)i{x) : 

^i{k) = —= / e-''"'Mx)dx. (1.3) 



27r 

Then the (unique) solution of the real elliptic Dirichlet boundary- value problem (1.1,2) 
above is given by 

In common with any solution of an elliptic partial-differential equation with analytic co- 
efficients, (/) is automatically a (real- or complex-) analytic function of both arguments r 
and X [9]. 

One can then rotate the Euclidean proper-distance-at-infinity into the complex: T — > 
rexp(i6'), where < 9 < it/2. The Laplace equation (1.1) still holds, although it is 
perhaps more natural to view it with respect to coordinates (ye^x), where yg — Texp{i9), 
in which case the potential is, for each fixed 6, a (complex) solution of the complexified 
Laplace equation 

e«I^.S = 0. (1.5) 



dye^ dx'^ 

Provided that < 9 < n/2 , this differential equation is strongly elliptic in the language of 
[11]. We continue to pose the same boundary data (1.2), with cf) — on the initial surface 
yo = , and with 

(f>{y0 = T,x) = 01 (a;) (1.6) 
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on the final surface yg — T [equivalently given by r = Texp(— i^)]. The corresponding 
complexified solution (f){yg,x) of this boundary-value problem still exists and is unique, 
provided always that < 9 < it/2 , being given by the analytic continuation of Eq.(1.4): 




= ye exp 



\ 1 f°° 



ikx 



sinh(A; J/51 exp(— z^)) 



$i(A;) dk . 



(1.7) 



sinh(/cTexp(-z^)) 



Provided that < 9 < tt/2 , the denominator in the integrand of Eq.(1.7), namely 

sinh(A;Te-*^) = cos (/cT sin ^) sinh(A;T cos ^) - z sin (/cT sin ^) cosh (/cT cos ^) , (1.8) 

is non-zero for all k^O . By examining a small neighbourhood of = , one sees that 
both numerator and denominator have simple zeros at /c = , whence the integrand in 
Eq.(1.7) is well-behaved for all k. In this way, it can be shown that the integral (1.7) 
provides the (unique) solution of the boundary- value problem, for < 9 < 7r/2 . 

On the other hand, for the exactly Lorentzian case with 9 = 7r/2 , one can see that 
no solution exists for the Dirichlet boundary-value problem, for typical choices of the final 
boundary data (t>i{x), or equivalently of its Fourier transform $i(/c). This is because, 
following the representation (1.4), one now has sin(/cT) for the denominator, rather than 
sinh(A;T). The existence of a (smooth) Lorentzian solution (p{t, x) would then imply that 
$i(/c) vanishes for all those values of k , namely, for all A; = nyr/T (n a positive or negative 
integer, excluding the case n = 0), such that sin(A;T) vanishes. A typical </>i(a;) will not 
have $i(?i7r/T) = for a single value of n (integer, ^0). This argument can be refined, to 
show rigorously the non-existence of a solution to the above (non-trivial) Dirichlet problem 
for the fiat wave equation. Thus, in this example, the Dirichlet boundary-value problem 
for the wave equation is badly posed. 

In [1] and in this paper, we study the field-theoretic analogue of the boundary-value 
problem above, with data specified on an initial and a final spacelike hypersurface, with 
either a (Riemannian) proper distance at spatial infinity, T, between the surfaces, or a 
complex-rotated separation Texp(i^), where < < 7r/2 . For simplicity, we consider Ein- 
stein gravity coupled to a massless scalar field (f) , and first restrict attention to spherically- 
symmetric configurations, as giving a large number of 'background' or 'reference' space- 
times describing gravitational collapse to a Schwarzschild black hole in Lorentzian signature 
[12,13]. The metric is taken (without loss) in the form (2.1) below; the corresponding scalar 
field is of the form 0(t, r). In Riemannian signature, the Einstein field equations are given 



As yet, very little is known rigorously in the theory of partial differential equations 
about existence and uniqueness for this boundary- value problem - even for the purely real 
Riemannian case with a spherically-symmetric 3-metric (/iy)/,F and scalar field {(p)i,F 
specified on the initial and final hyper surfaces, and a real proper-distance separation T 
given at spatial infinity. In the limit of weak boundary data, and hence, via the classical 
Riemannian field equations, also of weak fields in the interior, one would expect fixed- 
point methods, analogous to those used in [15], to establish existence and uniqueness for 
the Riemannian boundary-value problem. This weak-field Riemannian problem has also 



in [14]. 
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been investigated numerically [14] . For weak scalar boundary data, global quantities such 
as the mass M and Euclidean action / appear to scale quadratically, in accordance with 
analytic weak-field estimates [16]. 

Considering still the real Riemannian case, there does not 'experimentally' appear to 
be any obstacle to extending the above (good) results numerically to strong boundary data, 
except for the ever-larger amounts of computer time required. It may be that a typical 
pattern will emerge numerically for the general 'shape' of the classical gravitational and 
scalar fields, in the limit of strong-field boundary data. In that case, it might be possible 
to find analytic approximations for the strong-field limit (quite different from those valid 
in the weak-field case), which could provide a large amount of further analytical insight 
into the solutions of the strong-field Riemannian boundary-value problem. It is, of course, 
such strong-field boundary- value problems which are of most interest, when we consider 
rotating: T — > Texp(z^), with < ^ < 7r/2 , towards a Lorentzian time-interval. 

The strong-field limit for this boundary-value problem is expected to correspond, in 

the Lorentzian limit in which 6 approaches 7r/2 from below (although 6 never attains 
the value 7r/2), to classical spherically-symmetric Einstein/scalar solutions which form a 
singularity, surrounded by a black hole. As described in the Abstract above, one expects 
that arbitrary smooth gravitational and scalar data can be posed on the final spacelike 
hypersurface E p , allowing for a rotation into the complex of the Riemannian distance- 
separation T (> 0) between E/ and Hp , oi the form T — * Texp{i6), with < 9 < tt /2 , and 
provided that "Ep has the correct mass M. The (real) magnitude T could then be chosen 
arbitrarily large, corresponding to a non-singular complex solution of the Einstein/scalar 
field equations between E/ and a final hypersurface Ej? at late time. In this sense, in 
the complexified solution, one expects to have Bogoliubov transformations to surfaces "Ep 
which are 'after the singularity in the Lorentzian-signature geometry', provided always that 
one does not take the limit 6 = 7v/2 . In the same sense, one may say that 'the singularity 
is simply by-passed in the analytic continuation'; by analogy with the scalar- field example 
of Eqs. (1.1-8) above, there should be an analytic complexified classical solution to the 
boundary- value problem, which only reaches a singular boundary precisely at Lorentzian 
signature [9 = 7r/2). Provided that one retains the freedom to rotate T into the complex 
(as was insisted upon also by Feynman [8] in his +ie prescription), one has the possibility 
of circumventing Lorentzian singularities by deforming time-intervals suitably into the 
complex, much as one avoids singularities of functions f{z) in the ordinary complex z- 
plane. 

For the more general Einstein/scalar boundary- value problem, in which one allows for 
weak non-spherical perturbations in the boundary data (hij, 4>)i^f , the (complex) classical 

action S'dass can be studied as a functional of the perturbative final data (/ij-j'', 4>^^'')f and 
a function of the complex variable T, if for simplicity we regard the spherically-symmetric 
initial data (/ly, as fixed. The classical action .Sdass can be calculated, and thence the 
semi-classical amplitude, proportional to exp(i S'dass)- Finally, the limit ^ ^ 0+ can be 
taken, yielding the Lorentzian amplitude. 

This description is clearly very different from that usually adopted for black-hole 
evaporation [17-20], which involves the study of Bogoliubov transformations between an 



4 



early- and a late-time surface, where the late-time surface passes inside the black-hole event 
horizon, close to the curvature singularity. The standard approach yields probabilities for 
configurations of outgoing particles, and a late-time density matrix, but not quantum 
amplitudes for states defined typically by weak data over an M"^ surface at late times, as 
in the present case. It is intended, in this Letter, to make some contact between these 
two approaches, even though our approach addresses a wider class of questions than does 
the usual approach. In particular, we show how the Bogoliubov transformation between 
our two R'^-surfaces E/ and "Ep can be described in terms of the quantities appearing 
in [2,3] (see also [1]). The Bogoliubov quantities \(3su!em\ for the original calculation of 
black- hole radiance [17] agree with those found by the above boundary- value approach, 
to high accuracy (Sec. 3). In Sec. 4, we consider the resulting probability distribution for 
configurations of the perturbative scalar field (p^^^ on the final surface Tip. It is found 
that this probability distribution cannot be expressed in terms of the diagonal elements of 
some density-matrix distribution; thus, the correspondence between our approach and the 
density-matrix approach of [18] is limited. 

2. Bogoliubov transformations between surfaces 

As in [1-4] , we write the Lorentzian-signature classical 'background' metric in the form 
ds^ = - e^^^'^-^de + e«(*'^)dr2 + r'^idO'^ + sui^edcj)^) . (2.1) 
For later reference, we also define the 'mass function' m{t, r) by 

2m(t, r) 



exp 



{-a(t,r))=l-'-^. (2.2) 



In the boundary-value problem outlined above, let us write {"^^v , $) for the 'background' 
spherically-symmetric metric and scalar field. Writing also V^^ for the background covari- 
ant derivative, consider a linearised classical solution 0^^^ of the scalar wave equation 

V^V^^^^) = , (2.3) 

Consider a basis {fu;'im{x)} of (nearly-) separated mode solutions adapted to S/ and 
another such basis {pojem{x)} adapted to Tip ■ On S/, the {foj'em} are chosen to be an 
orthonormal, complete family of complex solutions of the wave equations, which contain 
only positive frequencies {uj' > 0). Correspondingly for the {pcjim} on Tp. In a standard 
fashion [19], one expands out (l)^^\x) in terms of one or the other basis, as 

°° ^ /"OO 

0(1) (a;) = X] XI / ^<^'irn fuj'lm{x) + C.C.J , (2.4) 

°° ^ /"OO 

(/,(i)(a;) = ^ ^ duj b^im Pojimix) + c.c. . (2.5) 
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Note that we shall assume throughout that the asymptotic time interval T = \T\ ex.p{—i9) 
has negative imaginary part. For convenience of visualisation, one might wish to regard 
the angle 9 as being extremely small. Nevertheless, one only takes ^ ^ 0+ at the end of 
the calculation. Thus, this Section deals with a larger class of Bogoliubov transformations 
than those appearing in the usual formulation. 

The bases above may be chosen to have the form 

cl>u,lm{x) = N{u;) e-^-* YU^) , (2.6) 

r 

where N{u!) is a normalisation factor. On S/ , we require (f)ujem oc e~'^^'" at large r, where 
V = t + r* and where r* ~ r* , with 



* 

r 

' s 



r + 2M ln(^(r/2M) - (2.7) 



denoting the Regge- Wheeler radial coordinate [21,22] in an exactly Schwarzschild geometry 
of mass M. On , we require 0a;£m «: e~*'^" at large r, where u = t — r* . Thus, the 
{foj'em} are purely ingoing at infinity, while the {puiem} are purely outgoing. 

Since the {fuj'im} give a complete orthonormal set on S/ , a typical basis function 
(solution) Pu,£m on E/ may be expanded as 



Pojim — J dio' (oiuj'ujlm fui'im + Pu;'uj£m f^'£,-rnj ; 



(2.8) 



The sets of complex numbers {auj'ojem} and {Pu'wim} are the Bogoliubov coefficients, and 
the standard properties of Bogoliubov transformations are, for example, given in [13]. 

In our treatment [1-5] of linearised scalar perturbations (j)^^\x), propagating at late 
times through the Vaidya-like region containing outgoing black-hole radiation, almost all 
modes are adiabatic. That is, a general solution of Eq.(2.1) in this region can be built 
up from separated solutions of the form (2.6), with uj = —k and Ru,e{r) ~ Cke{f) • Here, 
Cki{r) obeys 



where, as usual. 



+ (A;' - V() = , (2.9) 



A = e(-«)/^^ , e-.e^.l-^I^, (2.10) 

or or r 



and 

Vdt.r) = - 

r 

The 'left' boundary condition on as r — > is 



= ^(^(^ + 1) + ^!^'). (2.11) 



CfeKO) = . (2.12) 
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The 'right' boundary condition, as r — > oo , is 

ikiir) ^ Zki exp(i/cr*) + zl^ exp(-i/cr*) , (2.13) 

where, for each {k,£), zm is a dimensionless complex constant. We then wrote out (f)^^'^ in 
terms of 'coordinates' {aktm} on the final surface Tip , 

1 oo £ »oo 

= - J2 J2 dk akem^keir)Yimm . (2.14) 

In the nearly-Lorentzian case, the late-time behaviour of quantum amplitudes [1-5] 
will be dominated by the eigen-frequencies 

kn = ^. (2.15) 

Correspondingly, one discretises the frequency integral (2.11). It will also be convenient 
to define real functions {fn£{f)} for n = 1, 2, . . . ; £ = 0, 1, 2, . . ., such that 

U{r) = {zn£ e'^-^' + 4, e-'^--") Uir) , (2.16) 

where /n^(r) ^1 as r ^ oo , fn£{r) ^0 as r ^ . 

It is now possible to make a comparison of this treatment of the quantum amplitude 
from initial data on to final data on M^, with the corresponding description in terms 
of Bogoliubov transformations (again, in our formulation, not in the standard formula- 
tion). One takes a discretised version of Eq.(2.5), appropriate to the final surface Ef , and 
compares the expansions of (j)^^^ in basis functions on Sf • In particular, one finds that 

\bn£mf = '^'^ kn\Zn£f \an£m + 0'-n£m\'^ = f£m{kn) , (2-17) 

where {bn£m} are the coefficients appearing in Eq.(2.5). We also find an expression for 
|6tj£rn| in terms of Bogoliubov coefficients. Including an implicit summation over the index 



m: 



^uj£m\ 



/■oo <>oo / 

2 J du' j duj" Re\/3^»^i al,,^g Cuj'£m Cuj"£-m 



(2.18) 



3. Particle emission rates for a final surface 



The total energy of the late-time Einstein/scalar system must be equal to M, the ADM 
mass [22] or energy of the original configuration. The 'second-variation' contribution E^'^^ 
to M, namely the contribution from all the particles produced, is 

/•oo 

E^^^ = Yl / d^^K£m\^ ■ (3.1) 
£m 



This can also be written in terms of the 'harmonic-osciUator coordinates' {cu,£m} of 
Eq.(2.4), on using Eq.(2.18). Further, the expected number < n^gm > of particles emitted 
[5] is 

< n^lm > = IPcoll , (3.2) 

2 

where \Puje\ is defined by 

/•CXD 

/ duj' Puj'uje P^'uj"£ = \Pu,ef S{u;,u") . (3.3) 
Jo 



In the original, standard, calculation of Bogoliubov coefficients and of probabilities for 
particle emission by a non- rotating black hole [17], neglecting back-reaction, it was found 
that ^ 

IPs^emf = rsu,im{m) (^e^^'^-i-lf'^ , (3.4) 

where rst^^^(m) is the transmission probability over the centrifugal barrier [23] and rh = 
2M(jj . This calculation, of course, referred to the case where one did not have a final 
surface Ei? of topology M^, with a freely-chosen metric on it, subject to the mass being M. 
But, because of the very-high-frequency (adiabatic) method in which this expression for 
\Psu)£m\ was calculated, it should still be valid (up to tiny corrections) in our present 
case, noting that the calculation of scattering over the centrifugal barrier should allow for 
the extremely slow time- and radial-dependence of the mass function m{t, r) in Eq.(2.9-ll). 



4. Probabilities and density matrices 



Following from the Lorentzian quantum amplitude calculated in [1-5] , for perturbative 
spin-0 data given on the final surface Ei? , while taking (for simplicity) exactly spherically- 
symmetric 'background' data on E/ , one finds the corresponding conditional probability 
density over the final perturbative scalar boundary data (j^^^^Yip ' given that (j^^^^lj^^ = . 
For a large (real) Lorentzian time-interval T, this is 



{O'kem}', \T\ 



= Ne 



-5\T\Mi 



exp 



{C'kem}', \T\ 



(4.1) 



is calculated 



where is a suitable normalisation factor. The quantity 'S'^ia.gg {o-kim}', T 
in [3]. 

The probability distribution (4.1) arises, of course, from squaring a quantum ampli- 
tude. But, given the thermal nature of black-hole evaporation in the usual description, one 
can still ask whether the probability distribution can be viewed in terms of the diagonal 
elements of some density matrix in a suitable basis [18-20]. This is simpler to understand 
if one assumes the relations 



(4.2) 
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/•OO 

^0 



(4.3) 



These conditions (4.2,3) do hold for the steady-state Bogohubov coefficients in the calcu- 
lation which neglects back-reaction on the metric [17]. One needs to check whether this 
diagonal form persists under adiabatic propagation through a slowly-varying potential. 
The slow variation of the potential with t and r occurs, as above, because, in the Vaidya 
metric [6] describing the radiative part of the space-time, the Schwarzschild mass M is 
replaced by the mass function m{t,r) of Eq.(2.2,10). 

We shall need the relation [19] 



(4.4) 



From these equations, one finds 



{(^uiim} 



= P 



{Cu>£m} 



= iV exp 



nlm 



' 1 + 2 \(3nl\^ ] \Cnlm'^ 



(4.5) 



As expected in the linearised theory, the modes evolve independently, corresponding to the 
product over nlm. 

In studying the question of whether P {c^£m} can be related to a density- matrix 
distribution, we shall need the relation 



exp 



OO 

-a;s/(l-s)] = ^Lk{x)s'', |s| < 1 , 



k=0 



(l-s) 

which gives the generating function for the Laguerre polynomials {Lk{x)} [24]. We set 



(4.6) 



< 1 



and 

X = 2 Xn£m = 4: (AcOn) \Cnem\ , 

which is dimensionless. Writing j = n£m, we find 



(4.7) 



(4.8) 



n 


exp(- 


-Xj) exp(-2|ft|^X,-) 


n 


exp(- 


OO 

-X,) Y,P{k,)L,^{2X,) 






kj=0 


E 

kj 




) n exp{-X,) L,^{2X,) , 

3 



(4.9) 



where P{kj) is defined as 



1 



(3, 



2k. 



nim 



Pikj) - --2 - > (4-10) 

\aj\ Oij 

giving the (normahsed) probabihty to observe the field in the state kj. Equivalently, P{kj) 
is the probability of finding kj particles outgoing at future null infinity, in the mode (nim). 
Eq.(4.9) describes the probability distribution P[{cj}] for the final scalar field, with the 
help of the Bose-Einstein distribution (4.10). From Eq.(4.5), -P[{cj}] depends only on the 
modulus \cj\ of the complex number Cj ; hence, the sum in Eq.(4.9) is diagonal. 

In Eq.(4.9), the Laguerre polynomial Lk{x) has k real, distinct roots. For positive 
argument Xj > and excited states {kj > 0), the function exp{—Xj) Li^.{2Xj) takes 
negative values for certain ranges of Xj . Therefore, this function cannot be interpreted as 
a probability density. 

For comparison, consider a density matrix of the form 

p{x,y) = Wiipi{x) i;*{y) , (4.11) 

i 

where the {ilji{x)} form a complete set and Wi is the probability to be in the state '^i{x). 
If the system is in the single state defined by ipi{x), then the probability density for 
observing x is |'0j(x)|^ . There is not a probabilistic interpretation for the full p{x,y) , but 
the diagonal contribution 

P{x) = p{x,x) ^ Yl""^ (4-12) 

i 

is the probability density for observing the coordinate x . In our example, if the field is in 
the kj-th. state, then N ex.p{—Xj) Lk.{2Xj) is the 'probability density' for the 'harmonic 
oscillator coordinate' Xj . We then obtain the statistical 'probability' distribution (4.9) 
on multiplying by the probability P{kj) to be in the kj-i\v state, then summing over all 
possible states kj. This gives an interesting way of regarding the probability distribution 
(4.9), but it does not provide a sensible density matrix, based on the quantum amplitude 

exp S'^j^j^^gg [{afc^mjjT"] j of [3-5]. This shows that the density- matrix approach, based 
conceptually on linearised quantum field theory around a classical space-time, does not 
apparently replicate the results of a fully quantum approach (at least, not in the locally- 
super symmetric case). 



5. Conclusion 



Here, we have made some connection between the boundary-value approach of [1-5] 
to calculating quantum amplitudes for generic late-time configurations defined (say) on 
M"^, which involves a rotation of the time-at-infinity: T \T\ ex.p{—i9) into the complex, 
and the familiar Bogoliubov-coefficient approach to calculating probabilities in black-hole 
evaporation, with a final hypersurface near the space-time singularity. In the case that 
both the initial and the final space- like hypersurfaces E/, Si? have the topology of M^, there 
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is a relation between the Bogoliubov transformations in the two different approaches. The 
|/3st^£m| quantities for the original (Bogoliubov) calculation of black-hole radiance and 
for the topologically-M*^ approach agree, at least very closely. However, the probability 
distribution calculated from the 'M"^' quantum amplitude cannot be expressed in terms of 
the diagonal elements of some density-matrix distribution, so supporting the view that the 
two approaches are genuinely different. 
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